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IMPLEMENTING  BRENT’S  MULTIPLE-PRECISION  PACKAGE  ON  THE 
TEXAS  INSTRUMENTS  ADVANCED  SCIENTIFIC  COMPUTER 

USING  BRENT'S  MULTIPLE- PRECIS ION  PACKAGE  ON  THE  ASC 

This  note  will  enable  an  ASC  user  to  easily  and  quickly  use  MP,  a 
multiple-precision  arithmetic  package  written  by  Richard  P.  Brent  of  the 
Australian  National  University  [l].  This  note  is  intended  as  an  introduc¬ 
tory  guide  for  using  MP  on  the  Texas  Instruments  (TI)  ASC,  and  is  not  meant 
to  replace  the  documentation  for  MP,  Documentation  for  MP  is  catalogued  at 
USERCAT/ LIB/ UNSUP/ EXTDOC;  additional  documentation  and  source  listing  reside 
at  USERCAT/ LIB/UNSUP/EXTSCR.  The  object  library  (FX  compilation)  for  MP  is 
catalogued  at  USERCAT/ LIB/UNSUP/ EXTOBJ  and  contains  100  subroutines.  Refer¬ 
ence  [l]  is  recommended  for  those  who  would  like  to  make  a  deeper  study  of 
the  mathematics  of  the  MP  package. 

The  following  lists  the  minimum  JSL  for  a  typical  batch  job: 

/  JOB.... 

/  ASG  OBJ,  USERCAT/ LIB/ UNSUP/ EXTOBJ,  USE  =  SHR 
/  FTN 

fortran  (main)  program  and  subroutines  (if  any)  » 

/  LNK  LSPACE  =  XXXXX 
LIBRARY  OBJ 

/  FXQT  j 

/  EOJ 

The  following  non-executable  statements  must  be  included  in  the  calling 
program: 

COMMON  B,  T,  M,  LUN,  MXR,  R 


INTEGER  B,  T,  R 

R  must  be  dimensioned  and  B,  T,  M,  LUN,  MXR  must  be  defined. 

i.  For  best  efficiency  and  space  utilization  set  B=16384. 

ii.  Error  messages  are  printed  on  logical  unit  LUN,  where  LUN  is 
set  by  the  user;  normally  LUN=6. 

iii.  If  N  is  the  "equivalent"  number  of  floating  decimal  places 
desired  by  the  user,  then  T  is  the  smallest  integer  such  that 


For  B=16384, 


T  >  max  |  1+7  (log^B)"1,  l+N(log1QB)" 1  j. 

T  >  max  |  3,1+  .24N  j. 
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iv.  Exponents  can  lie  in  the  closed  interval  [-M, M]  ,  where  M  is  set 
by  the  user.  M  must  be  greater  than  T  and  less  than  or  equal  to  the  integer 
part  of  0.25  (231  -  1)  =  536,870,911. 

v.  The  value  of  MXR  is  obtained  as  follows:  For  each  subroutine  in 
the  MP  package  is  associated  a  parameter  called  SPACE,  which  is  a  linear 
function  of  T;  i.e.,  SPACE  has  the  form  SPACE  =  c*T  +  8 ,  where  a  and  0  are 
integers.  Suppose  that  n  subroutines  belonging  to  the  MP  package  are  used. 
Then,  for  each  of  these  routines  we  have: 


SPACE . 
J 


V  + 


1  £  j  S  n. 


Suppose  that  the  value  of  integer  T  as  defined  above  has  been 

chosen,  i.e.,  T  =  T  .  Then  MXR  is  defined  by 
o 


max 
1=£  j^n 


SPACE . 
J 


1  = 


max 

l^j^n 


a .  T  +  3  . 
J  °  J 


vi.  R  must  be  dimensioned  in  the  calling  program;  its  dimension  is 
MXR.  This  array  is  used  as  work  space  for  the  collection  of  MP  routines 
used. 


We  note  that  the  above  parameters  in  common  may  be  varied  dynamic¬ 
ally  subject  to  the  restrictions  of  the  MP  package  and  ASC  FORTRAN.  Further, 
we  note  that  MP  uses  "T-digit  floating  point  numbers"  which  are  stored  in 
integer  arrays  of  dimension  T+2  (see  documentation  and  example  programs). 

It  is  important  to  note  that  the  subroutine  MPSET  must  not  be  used 
as  it  will  produce  integer  overflow.* 

TABLE  OF  BASIC  ARITHMETIC  ROUTINES  AND  SAMPLE  PROGRAMS 

There  are  100  routines  in  the  MP  package.  Some  basic  routines  are 
[l,  pp  70  -  74  or  catalogued  documentation  previously  mentioned]: 

MPADD  Adds  two  MP  numbers,  SPACE  =  T  +  4 

CALL  MPADD  (X,  Y,  Z)  means  Z  =  X  +  Y 


MPSET  causes  integer  overflow  when  computer  wordlength  is  less 
than  48  bits.  Brent  suggests  two  alternatives  for  MPSET;  a  rewriting  of 
this  software  or  setting  the  parameters  B,  T,  M,  LUN  and  MXR  without 
calling  MPSET.  The  latter  has  been  chosen  since  it  provides  greater 
flexibility.  MPSET  was  not  removed  from  the  MP  package  because  it  is 
mentioned  (not  used)  elsewhere  in  the  documentation. 
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MPSUB 


Subtracts  one  MP  number  from  another,  SPACE  =  T  +  4 


CALL  MPSUB  (X,  Y,  Z)  means  Z  =  X  -  Y 
MPMUL  Multiplies  two  MP  numbers,  SPACE  =  T  +  4 
CALL  MPMUL  (X,  Y,  Z)  means  Z  =  X  *  Y 
MPDIV  Divides  two  MP  numbers,  SPACE  =  4T  +  10 
CALL  MPDIV  (X,  Y,  Z)  means  Z  =  X/Y 

MPIN  Converts  fixed-point  number  read  under  A1  format  to  multiple 
precision,  SPACE  =  3T  +  11 

MPOUT  Converts  multiple-precision  to  a  form  suitable  for  printing  under 
A1  format  (corresponds  to  F  or  1  formats),  SPACE  =  3T  +  11 

Appendix  1  gives  a  detailed  program  illustrating  how  these  arith¬ 
metic  routines  may  be  used  on  the  ASC.  Program  output  also  provided. 

An  example  program  TSTMP  using  MP  subroutines  is  also  given  in 
Appendix  2  of  this  note;  the  output  is  provided.  This  example  program  is 
an  ASC  adaptation  and  extension  of  the  example  program  supplied  with  the 
MP  package.  The  program  TSTMP  uses  (implicitly  and  explicitly)  49  of  the 
99  usable  MP  subroutines.  The  times  for  computations  are  included. 

MATHEMATICAL  FUNCTIONS,  CONSTANTS  AND  SOFTWARE 

The  capabilities  of  MP  include  the  following  [l,  p.  58-59]: 

i.  conversion  of  integer,  real,  and  double-precision  numbers  to 
multiple-precision  format  and  vice  versa; 

ii.  multiplication  and  division  of  multiple-precision  (mp)  numbers 
by  small  integers; 

iii.  addition,  subtraction,  multiplication,  and  division  of  mp  num¬ 
bers  ; 

iv.  powers  and  roots  of  mp  numbers; 

v.  elementary  functions  of  mp  numbers:  arcsin,  arctan,  cos,  cosh, 
exp.  In,  sin,  sinh,  tan,  tanh; 

vi.  some  special  functions:  Bessel  functions  of  the  first  kind, 
error  and  complementary  error  functions,  exponential  and  logarithmic 
integrals,  Dawson's  integral,  and  gamma  function; 

vii.  constants:  Bernoulli  numbers,  rr,  Euler's  constant,  and  the 
Riemann  zeta  function  for  positive  integer  arguments; 

viii.  fixed  and  floating-point  decimal  output  and  free- field  decimal 
input  of  mp  numbers; 


i 
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ix.  integer  and  fractional  parts  of  mp  numbers; 

x.  routines  for  error  handling,  testing,  and  debugging; 

xi.  miscellaneous:  comparison  of  mp  numbers,  storing,  packing  and 
unpacking  mp  numbers,  etc, 

COMPUTATION  TIME  COMPARISONS 

Computation  times  for  the  basic  arithmetic  operations  of  the  MP 
package  and  for  the  TI  package  [2]  are  given  in  Tables  1  and  2.  Table  3 
clearly  shows  the  superiority  of  MP  over  the  Texas  Instruments  Extended 
Precision  Package  [ 2 J  with  respect  to  each  of  the  arithmetic  operations. 

IMPROVING  MP  ON  THE  ASC 

Although  the  speed  of  the  MP  package  is  substantial  (see  timings 
listed  in  example  output),  the  speed  of  the  package  could  be  increased  by 
optimizing,  for  the  ASC,  the  basic  addition,  multiplication,  and  division 
routines : 


MPNZR,  MPMLP,  MPDIVI ,  MPADD2 ,  MPADD3 ,  MFMUL2. 
TABLE  1 


33  significant 

figures  + 

- 

X 

MP  package 

3,910 

4,090 

9,210 

77,430 

TI  package 

28,600 

42,100 

1,370,000 

290,000 

The  entries  of  this  table  give 
for  operations  on  two  multiple- 
figures.  TI  package  refers  to 

“8 

the  time  in  clocks  (1  clock  =  8x10  second) 

precision  numbers  with  33  significant 
reference  [2]. 

TABLE  2 

70  significant 

figures  + 

- 

X 

■+ 

MP  package 

5,570 

5,780 

14,200 

109,000 

TI  package 

52,400 

81,000 

5,820,000 

873,000 

The  entries  of  this  table  give  the  time  in  clocks  for  operations  on  two 
multiple-precision  numbers  with  70  significant  figures. 
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TABLE  3 


r 

+ 

- 

X 

+ 

33  significant  figures 

m 

10.3 

149 

3.75 

70  significant  figures 

9.4 

14 

410 

8 

The  entries  of  this  table  give  the  ratio  r  of  TI  Extended  Precision  Package 
to  MP  Package  for  each  operation  performed  on  ASC. 

CONCLUSION 

Brent's  multiple-precision  arithmetic  package  is  now  available  on  the 
ASC.  It  not  only  performs  the  basic  arithmetic  operations  up  to  orders  of 
magnitude  faster  than  what  is  now  available  on  the  ASC,  but  it  is  also 
supported  by  we  11- documented  software  of  functions  and  constants  of  great 
interest  in  scientific  computations. 
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